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INTRODUCTION 


The design of future large space structures will require the development 
of advanced control systems in order to meet stringent pointing requirements. 
There are two types of operational environments that are of interest. First, 
earth pointing applications from geostationary orbits. These applications are 
characterized by vibration and shape control problems which must operate within 
a persistent disturbance environment. Second, maneuvering or slewing 
applications where the vehicle is transferred from an initial pointing direction 
to a final desired pointing direction in finite time. A key issue for maneuver 
problems is that the desired rigid body pointing direction is achieved while 
suppressing any maneuver induced structural response. The technical solutions 
for these applications are fundamentally different. For example, the vibration 
and shape control applications represent infinite time or steady state problems, 
where the major concern is system stability. The stability questions arise 
because the control system is continuously generating command inputs in an effort 
to suppress the structural response induced by the persistent disturbance 
sources. As a result, if the mathematical model used for designing the control 
has errors, the potential exists for exciting a resonance response due to control 
spillover effects. Slewing control problems, on the other hand, represent 
finite-time control problems because the maneuvers are usually designed to be 
accomplished within a specified period of time. Even though the maneuver loads 
tend to be large compared with typical on-board or environmental disturbances, 
their finite time duration tends to minimize the potential for exciting resonance 
responses due to modeling errors. 

The report presents an advanced control design problem and discusses the 
requirements for automating both the structures and controls design efforts. 
First, a general three dimensional slewing technique is developed and applied to 
a large geostationary platform. This system presents significant technical 
challenges because of two large flexible antennas which are attached to a 
flexible truss structure. The control problem formulation treats the full 
nonlinear coupled rigid and flexible body control problem. Second, the software 
development issues are considered for developing an integrated structures and 
control design environment. This development is driven by the recognized need 
to coordinate both the structures and control design efforts to minimize mission 
threatening control/structure interactions. 
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SLEW CONTROL METHODOLOGY 


Control Strategy 

The structure used in this study is the Earth Pointing Satellite (EPS) 
shown in Fig. 1. The control methodology for the Class III slew control of the 
EPS model follows the approach of Chun, Turner, and Juang, 1987. The approach 
consists of an optimal open-loop maneuver designed for a rigid model of the EPS 
structure, and a closed-loop perturbation feedback designed for a reduced order 
flexible model. The rigid body controls are applied to the flexible plant. The 
difference between the nominal rigid body response and the actual flexible 
response is then fed to the perturbation feedback controller to generate control 
corrections that cause the flexible structure to follow the rigid body response 
as closely as possible. Both the open-loop and the closed-loop controllers are 
designed using frequency- shaping techniques that impose a two-pole roll-off in 
the frequency domain of the control effort. 

While the results of Chun, Turner, and Juang, 1987 were obtained with a 
fair separation of control and structure bandwidths, the current study involves 
maneuvers where the bandwidths are much closer. The current work also introduces 
a disturbance accommodating feedforward term into the control strategy that 
compensates for the excitation of the elastic degrees of freedom due to the 
nominal controls. 

Model Development 

The structure is developed by the CSI Analytical Design Methods (ADM) team 
at NASA LaRC. The structural design is a derivative of a geostationary platform 
designed in support of Earth Observation Sciences (EOS). 

The total mass of the vehicle is 1028 kg. It consists of 548 kg of 
structural mass, 329.63 kg of nonstructural mass on the antennas, and 150 kg of 
actuator mass. The truss in Fig. 1 is about 25 m long while the two antennas are 
15 m and 7.5 m in diameter. There are supports for each of the antennas. A more 
detailed description of the structure and material properties are presented in 
the following paragraphs. 

The truss is a 3 meter total cross section structure. The truss consists 
of 51 mm diameter x 1.59 mm wall thickness graphite tubes. The modules of 
elasticity and mass density for the tubes are 275.9E+09 N/m 2 and 3250 kg/m 3 
respectively. There are total of 135 beam elements that form the truss. 
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Both the 15 m and 7.5 m antennas have 12 ribs which radiate outward from 
the center. The ends of the ribs are connected by hoop elements. Graphite tube 
diameters and thicknesses for the antennas are the same as those used for the 
truss elements. The support of each antenna consists of four beam elements that 
offset the antenna from the main truss (see Fig. 1). 

Torquers are located at the CG of the vehicle, acting about all three axes 
to provide the slewing motion. The support for the actuators is made of the same 
element sizes and mechanical properties as the truss elements. 

The first 20 undamped vibration frequencies and their corresponding modal 
descriptions are: 

Mode No. Frequency in Hz Description 


1-6 

0.0 

7 

0.2422 

8 

0.4056 

9 

0.5652 

10 

0.6562 

11 

0.8882 

12 

0.8882 

13 

1.4377 

14 

1.5360 

15 

1.7762 

16 

1.7762 

17 

3.0258 

18 

3.0258 

19 

3.1532 

20 

3.1532 


Rigid Body 

Torsion of 15 m Antenna 
z- rotation 15 m Antenna 
x- rotation 15 m Antenna 
Torsion 7.5m Antenna 
15 m Antenna mode 
15 m Antenna mode 
z-rotation 7.5 m Antenna 
x- rotation 7.5 m Antenna 
15 m Antenna mode 
15 m Antenna mode 
15 m Antenna mode 
15 m Antenna mode 
15 m Antenna mode 
15 m Antenna mode 


It should be noted that the first 20 undamped elastic modes are all antenna 
modes. The first truss vibration mode occurs at a relatively high frequency, 
about 7.0 Hz, which is the 50th elastic mode. The NASTRAN EPS model has been 
generated by Joe Walz at NASA LaRC. It consists of a total of 95 nodes and 223 
elements. The truss and the two antennas are treated in the model as a single 
body. The finite element data used in DISCOS includes the modal frequencies of 
the first 45 elastic modes, the translation and rotation modal shapes for all 95 
nodes, and lumped masses at each node point. Modal damping of 0.5% is added for 
each mode. Inertial integrals are generated by DISCOS from mode shape and lumped 
mass data to account for the coupling between rigid and flexible body dynamics. 
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Actuator Placement 


It is well known that force actuation at various locations of a large space 
structure can provide significant vibration suppression of the structure. A 
preliminary force actuator analysis is performed to determine the force actuator 
locations. The objective of the actuator placement task is to find a set of 
feasible, but not necessarily optimal, actuator locations so that force actuators 
can be used for performance enhancement. 

Potential force actuators can be lightweight variable cold gas thrusters, 
and potential torquers can be CMG or reaction wheels. The force actuators and 
torquers are modeled as external forces and torques; no force actuator, CMG or 
reaction wheel dynamics are modeled. The commanded forces and torques are 
frequency- shaped so that there is very little energy at the high frequency end 
of the control spectrum. Furthermore, the force actuators and torquers are not 
output limited. 

For this study, two actuator placement algorithms have been examined - an 
H 2 optimization algorithm (Fanson, Blackwood, and Chu, 1989) and a degree of 
controllability algorithm (Longman and Alfriend, 1990, and Viswanathan, Longman, 
and Likins, 1984). The H 2 optimization algorithm was developed by Jim Fanson at 
the Jet Propulsion Laboratory to preselect the actuator locations on a truss 
structure. This algorithm has the advantage that it is theoretically simple and 
numerical computation can be readily performed by standard control design tools 
such as Ctrl-C. Its disadvantage is that the algorithm lacks strong theoretical 
support. The H 2 formulation assumes that the external disturbances are white or 
colored noise. 

Based on engineering intuition, one would expect elastic deformations at 
extremal locations of a large spacecraft to be large. Hence, the outer rims of 
the EOS antennas are expected to provide better actuator locations than locations 
on the truss. The suitability of the locations at the outer rims of the antennas 
has been confirmed by both simulation and actuator placement studies. A detailed 
description of the algorithm is presented in the following paragraphs. 

The cost function is selected to be a linear combination of the total 
energy of the position vector y at specified nodes and the control vector u, 
subject to disturbances generated by the torquers located at the CG. It is 
assumed that the disturbance is a low bandwidth colored noise. This formulation 
provides a well-posed H 2 optimization problem where an analytical closed- form 
solution exists. Full state information is assumed in the present analysis. 
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In mathematical terms, the H 2 optimization problem for actuator placement 
is defined as 


P - min II z II? 
VMO,®) 


( 1 ) 


where 


z = [77 u] T ; 7>0 

7 is a positive constant weighting factor, y(t) and u(t) are the output and 
control of the system, respectively, and ||•|| represents the Euclidean norm. 

For simplicity, only the modal state equations are considered. The 
corresponding state- space equations are 


x = Ax + Bu + b M w 


( 2 ) 


and 



The unique optimal control is expressed as 

u = Fx 


(3) 


( 4 ) 


where 

F = -B 7 X 

and X is the positive definite solution to the following Riccati equation: 

A t X + XA - XBB r X + 7 2 C 7 C = 0 
The minimal cost can be expressed as 
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(5) 



The total energy in y and u can be computed by using 


j|y|| 2 = trace(cw c C T ) ^ 


( 6 ) 



trace(FW c F T )J 


V2 


where W c is the solution of the following Lyapunov equation: 


(7) 


(A + BF) W c + W c (A + BF) t + BB t = 0 


The total energy in each control input, uj, of the optimal control u is computed as follows: 


lhl 2 = (b i T W 0 b i ) V2 


( 8 ) 


where bj is the i-th column of the control influence matrix B, and W 0 satisfies the following 
Lyapunov equation: 


(A + BF) t W 0 + W 0 (A + BF) + F t F = 0 


In this preliminary actuator placement analysis, the disturbance input is provided by the 
torquers located at the CG (node 95) of the vehicle. The candidate force actuator locations are 
arbitrarily selected to be nodes 69, 66, 63, 60, 52, 95, 94, 88, and 78 along all 3 axes, while the 
output of the system is chosen to be at nodes 69, 66 and 95 (CG of the vehicle) along all 3 axes 
(see Fig. 2). Equation (8) is used to evaluate actuator effectiveness. The five top ranking actuator 
locations based on a modal energy consideration are: 
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Rank 

Actuator 

Location 

1 

Node 

69 

X 

2 

Node 

66 

X 

3 

Node 

63 

X 

4 

Node 

60 

X 

5 

Node 

88 

X 


These five actuator locations are input into DISCOS for plant linearization 
as part of the perturbation feedback design. The linearized system dynamics 
matrix and the control influence matrix are used to compute feedback gains using 
LQR theory. However, the resulting closed-loop dynamics matrix is found to be 
unstable because the selected actuator locations do not provide the 
controllability required. A controllable system is obtained after switching two 
of these five actuators to be along the y and z directions, leading to: 


Actuator 

Direction 

Node 

69 

X 

Node 

63 

X 

Node 

66 

y 

Node 

60 

z 

Node 

88 

X 


Theoretically, the five force actuators, together with the torquers at the CG 
provide sufficient controllability. However, simulations indicate that the 
elastic displacements at opposite ends of the antenna rims tend to be opposite 
to each other in sign, with approximately the same magnitude. This observation 
suggests that the actuators be placed in pairs at opposite sides on the antenna 
rim where peak displacements are experienced. Accordingly, two force actuators 
are added to the above set of five actuators at node 66 in the z direction and 
node 60 in the y direction. No attempt has been made to add an additional 
actuator opposite to node 88 x because the magnitude of the elastic deformation 
at the 7.5 m antenna is about a factor of 5 smaller than that of the 15 m 
antenna. The final set of actuators selected is listed below: 


Actuator 

Direction 

Node 69 

X 

Node 63 

X 

Node 66 

y.z 

Node 60 

y.z 

Node 88 

X 


By including these force actuators, satisfactory performance is achieved. 
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I 



15 M Antenna 


Nodes 


Key results are presented in the Simulation Results section. 

A second actuator placement algorithm has also been studied; it is based 
on the concept of the "degree of controllability" (DOC). However, because of 
time constraints, no simulation results were obtained. 

Optimal Rigid-Body Open-Loop Control 

The rigid body equations of motion are given by a set of nonlinear 
differential equations, and the nominal rigid-body control is specified by the 
minimization of a finite-time performance index subject to specified initial and 
final conditions. The solution for the rigid-body open- loop control thus 
requires the solution to a nonlinear two-point boundary -value problem. The 
solution approach is identical to that described in Chun and Turner, 1988 and 
will only be briefly described here. 

While there are many methods for solving nonlinear two-point boundary-value 
problems (TPBVP) , the approach we have chosen employs the continuation or 
homotopy chain method, since the method is easy to apply,’ and it has global 
convergence properties. The continuation method involves the embedding of a 
parameter into the original problem to form a family of problems in such a way 
that a simple, easy- to-solve problem is defined when the parameter is set to 
zero, and the original hard-to-solve nonlinear problem is recovered when the 
parameter is set to one. In the version of the continuation method that we have 
chosen, the solution to the easy-to-solve problem is first obtained, and the 
continuation parameter is slowly incremented from zero to one, sweeping out a 
series of solutions which culminate in the desired solution of the original 
nonlinear problem. 

State Variables 


For the rigid body control problem, let us select angles, angular 
velocities, torques, and torque rates as the state variables. Here, we concern 
ourselves with only the rotational degrees-of-freedom, since translational motion 
is of little consequence for slew maneuvers. The available choices for angular 
representation are direction cosines, Euler parameters (quarternions) , or Euler 
angles. We choose Euler parameters because it is a unique set of numbers and 
does not encounter singularities at certain orientations, compared with a set of 
Euler angles, which is determined from one of twelve possible Euler sequences and 
exhibits "gimbal lock" at certain orientations. Direction cosines are also not 
chosen as state variables because of the larger number of variables needed. Body 
angular velocities are used as the velocity variables, and the torque variables 
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are normalized by the moment of inertia matrix. 


The Euler parameters (/9 0 , j9 lf /9 2 , 0 3 ) may be interpreted geometrically in 

terms of Euler's theorem which states: An arbitrary reorientation of a rigid body 
can be accomplished by a single rotation of the body about the principal line 2 
through the angle <f> . In terms of 2 and <f>, the Euler parameters are expressed as 


00 = cos -| 

= l i sin-^ , i=l,2, 3 


(9) 


It is clear from the above expressions that the Euler parameters satisfy the 
constraint 


E0i = i 


i-0 


( 10 ) 


and hence are a once -redundant set of variables. 


The Euler parameter rates are related to the body angular velocities by the 
equations 

j} = ( n ) 


where 


n(u) 



0 


-“V 


1 


0 


~ u y 

2 

u y 


0 

u x 


u z 

“V 

- w * 

0 


The differential equation for the body angular velocities is derived from 
Euler's equation, and is given by 

« = u 0 + G(u) (12) 


where 
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ff(u) = -[I] _1 [5][I]w 


and [I] is the moment of inertia matrix. 


The normalized torques and torque rates are included as state variables so 
as to incorporate frequency- shaping of the torque profile. The torque vector is 
normalized by the inertia matrix because this reduces the range of numerical 
values between the states and costates of the resulting optimal control problem. 
The normalized torque vector u 0 , which has units of rad/sec 2 , is related to the 
physical torque vector by the expression 


r -v 


[I] 


-il 


“x 

» 

“y 


where [u x , Uy, u 2 ] is the physical torque vector. The differential equations 
for the normalized torques and torque rates are given by: 

u 0 = u x . (13) 


and 



(14) 


Optimal Control Definition 


We define the optimal control problem for arbitrary frequency- shaped three- 
axis maneuvers of a general rigid body by the minimization of the following 
finite -time quadratic performance index: 



* ' 

J = |J,Nx T(C) “2^)^ 

x(t) 

u 2 (t) 

h 4 


(15) 


where the state vector is defined as 

X = [/3 T « T u* u’] T 

the control variable is defined as u 2 , and the weight matrix is defined as 
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"o 

0 



0 


0 0 0 0 
Q 0 0 0 

010 I/Uj 
0 0 0 0 
o i/4 o i/4 


The minimization is to be carried out subject to the state differential 
equations given by Eqs. (11)-(14), and specified initial and final states. The 
symbol I represents the (3x3) identity matrix. The form of the penalties on the 

normalized torque and torque rates reflects the time -domain equivalent of a 
frequency- domain penalty on the normalized torque (see Gupta, 1980). This 
penalty produces a two pole roll-off of the normalized torque's frequency content 
at the break frequency, Wg, and also allows the torque and torque rate to be 

specified at the initial and final times. A weighting penalty is placed on the 
angular velocities but not on the Euler parameters because the latter are allowed 
to be large. 

The costate differential equations can be shown to be given by: 


7 = 0(w)7 (16) 

X = -Qw - B t 7 + {[!«]} X [I] ml X - [I] [5] [I ] -1 A (1?) 

^0 = - u o ■ u z/4 - A (18) 

/*! “ -/*o (19) 


where 7(4x1), A(3xl), jj 0 (3x 1), and p x (3xl) are costate or adjoint variables for 
y9, w, u 0 , and u 1( respectively. The optimal control vector, which is the second 
time derivative of the normalized torques, is given by 

u 2 = + u 0 /wjj] ( 2 °) 
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Solution of the Nonlinear Optimal Control Problem 


The necessary conditions given by Eqs. (16) -(19) and the associated initial 
and final conditions represent a nonlinear TPBVP. The solution approach that we 
have taken employs the continuation method, whereby continuation parameters are 
embedded into the original nonlinear TPBVP to form a family of intermediate 
problems in such a way that a simple problem is defined when the parameters are 
set to zero and the original nonlinear TPBVP is recovered when the parameters are 
set to one. 


The simple problem that we have chosen is a single- axis maneuver about a 
principal axis of inertia, since its solution is known in closed form. To obtain 
this simple problem, we introduce a continuation parameter, a 1( that multiplies 


the off-diagonal terms of the inertia matrix: 

[!(«!>] = 


xx 1 xy 1 xz 




I yy ~ a i^yz 


- a l J zx -a i I zy 


which we will call the modified inertia matrix. We also introduce a second 
continuation parameter, a 2 , that multiplies the terminal angles, angular 

velocities, normalized torques, and torque rates along two of the three 
coordinate axes. The remaining coordinate axis is the selected axis about which 
the single-axis starting guess is to be computed. These will be called the 
modified boundary conditions. A modified TPBVP is defined when the modified 
inertia matrix is substituted into the nonlinear state-costate differential 
equations in conjunction with the modified boundary conditions. 


A separate set of minimum- state optimal control equations can be generated 
which correspond to the continuation-modified maneuver problem with the 
continuation parameters set to zero. The solution of the state-costate 
differential equations for this problem can be expressed in the exponential 
matrix form 


x(t t ) 




x(t D ) 


( 21 ) 


where x is the augmented state-costate vector and K is the Hamiltonian matrix. 
The linear TPBVP can be solved by the standard partitioned matrix solution, 
yielding the initial costates (Chun, 1986). The initial angle and initial 
angular costate can then be mapped into the corresponding two Euler parameters 
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and two Euler parameter costates for the continuation-modified problem with zero 
continuation parameters (p. 54, Chun and Turner, 1988). The other variables - 
angular velocities, normalized torque, normalized torque rate, and their costate 
variables are mapped directly. All other state and costate variables of the 
nonlinear problem are set to zero at the initial time. 

The continuation of the simple problem to the original problem then 
proceeds as follows: 

1. Increment and a 2 . Compute the corresponding modified 
inertia matrix and modified boundary conditions . 

2. Integrate the modified nonlinear state-costate differential 
equations and sensitivity matrices. 

3. Compute the terminal error. 

4. If the terminal error is large, apply corrections to the 
initial costates and go to Step 2. 

5. If o t - 1, and a 2 - 1, stop. The nonlinear solution has been 
found. 

6. Go to Step 1. 


Details on the calculation of the sensitivity matrices and the initial 
costate corrections can be found in Chun, 1986. 

Closed-Loon Perturbation Feedback Controller 

Since the open-loop optimal control is designed with a rigid model of the 
structure, direct application of the nominal control profile on the flexible 
structure will lead to errors in the maneuver angle as well as unwanted 
vibrations. To suppress these errors, a closed-loop perturbation feedback 

controller is designed to apply corrective torques and forces on the flexible 
structure to cause it to follow the nominal rigid body profile as closely as 
possible. The control design closely follows the approach described in Chun and 
Turner, 1988, with the exception of an added disturbance rejection to account for 
elastic response due to applied nominal torques. 

Control Design Strategy 

A reduced order model for the flexible plant is selected based on the 
response of the full order flexible plant to the open- loop nominal torques. 
Since the nominal torques are determined from a frequency- shaped control 
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formulation, the number of modes that are excited by the nominal torques is 
small, thus requiring only a small number of modes for the reduced order model. 

Since the plant dynamics is nonlinear during the large-angle rapid 
maneuver, we choose a piecewise-linear quasistatic approach to the control design 
so that standard linear quadratic regulator (LQR) techniques can be applied. The 
nonlinear flexible plant is linearized at several points along the nominal 
trajectory, and constant LQR feedback gains are calculated for each of the 
linearized plants. The number of linearization points to be used depends on the 
'strength' of the nonlinearities. 

Frequency- shaping is also used in the design of the feedback controller to 
minimize the possibility of control spillover. A disturbance- rejection 
feedforward is introduced to reduce the amount of excitation of the elastic modes 
that is caused by application of the nominal torques. Gain scheduling during the 
actual maneuver is performed by interpolating the LQR feedback and feedforward 
gains between the linearization times. 

Control Problem Formulation 


Flexible plant linearization is performed numerically by DISCOS (Dynamic 
Interaction Simulation of Controls and Structures), which uses finite differences 
and quadratic fits to obtain the linearized plant dynamics and control influence 
matrices. The linearization is performed about the nominal rigid body trajectory 
at several equally spaced intervals in time. 


For each linearized plant, a set of feedback and feedforward gains are 
obtained from a frequency- shaped LQR design. Specifically, we seek to minimize 
a performance index of the form: 


J = 



Sug(t) Su*(t) 


fSx(t) 1 


Su T 2 (t)] W 


Su 0 (c) 

' Su.Ct) ’ 

Su 2 (c) 

» - 


dc 


( 22 ) 


where 


W = 


Q N 
N J R 


Q 


Q 0 O' 
0 D 0 
0 0 0 


N J = [0 D / wj 0] Rm D/u* 
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Sx is a vector of rigid and elastic modal perturbations, fiu 0 is a vector 
containing actuator forces and perturbations of normalized torques, and S u lf Su 2 
are the first and second derivatives of Su 0 . The form of the weighting terms on 
the control-related states reflects the use of frequency-shaping which penalizes 
the frequency content of the controls that are above the break frequency, 

Minimization of the above performance index is carried out subject to the 
following differential equation constraints which describe the reduced order 
flexible structure and frequency -shaped controls: 

SX(t) =A6X(t) +B6u z (t) +d(t) (23) 


where 
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In the above expression, and B (i) are the linearized system dynamics and 
control influence matrices, respectively, and d(t) is a vector of arbitrary 
disturbances. 

Application of Pontryagin's Principle leads to necessary conditions 
consisting of the following state and costate differential equations: 


SX = [A-BR- 1 N t ]SX - BR‘ 1 B X X + d 


(24) 


X - -[A T -NR' 1 B r ]X - [Q-NR- 1 N t ]SX < 25 ) 

and the following optimal control equation: 


Su z - -R- 1 [B J X+N J 6X) 


(26) 
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One can assume a feedback/feedforward solution of the form 


A(t) = P(t)«(t) + £(t) (27) 

which, upon substitution into the necessary conditions, yields the following 
equations for the gain variables at steady state: 


P=0 = -PA-A T P + [PB*N]R' 1 [N t *B T P] - Q 


(28) 


£ = 0 = - [A-BR- 1 ^} T P+N T )] T e - Pd 


(29) 


The optimal feedback/feedforward control can then be given by the 
expression: 


fiu 2 (t) = -R- 1 (B r P sa + N T ')SX(t) 

+ R-'B 7 [A 7 ~(P SS B+N)R~ 1 B 7 ]~ 1 P ss d(t) 


(30) 


where the disturbance is assumed to be slowly -varying. 
Disturbance Vector 


The disturbance vector d(t) can be arbitrary. For the purpose of reducing 
the amount of slew- induced deformations, it is desirable to model the modal 
forcing term that is due to the nominal slew as a modal disturbance. From the 
equations of motion, it can be shown that the modal acceleration caused by the 
nominal slew torques is given by: 


a n - [o o i ] 1 
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(id 


(31) 


where the mass matrix has been partitioned into rotational, translational, and 
modal parts, and a is the matrix of modal slopes at the point of application of 

the torques, u. The disturbance vector d(t) thus contains one non-zero term, 
given by Eq. (31) in the appropriate row for modal disturbance. 
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Both 'constant' and time-varying disturbance vectors are tested in the 
example simulations. In the 'constant' case, the last term in Eq. (30) is 
evaluated at the linearization times and interpolated during the maneuver. In 
the time-varying case, the feedforward gains (coefficients of d(t) in Eq. (30)) 
are interpolated during the maneuver, but the disturbance vector itself is 
calculated at each moment using Eq. (31). 

Force/Torque Decoupling 

One of the observations that can made from the simulation runs was that the 
additional use of force actuators for the perturbation feedback controller leads 
to an increase in the torque levels and correspondingly large force levels . What 
appears to be happening is that the actuator forces produce moment couples which 
oppose the applied torques. One way of alleviating this is to utilize the force 
actuators to contribute to the required nominal torque profile (in addition to 
the torquers) using a torque-proportioning algorithm that reduces the amount of 
induced elastic deformation. An additional approach is to constrain the actuator 
forces so as not to produce any net torques. The control forces will then be 
able to affect only the flexible response (as well as translation) . This 
approach is briefly outlined as follows. 

For frequency- shaped controls, the decoupling constraint can be imposed on 
the highest derivative of the forces that appears in the control formulation 
(assuming zero initial conditions). The constraint to be applied can be 
expressed by the equation: 

E ? i f i ■ 0 (32) 

i 

where r A is the vector from the center of mass to the force actuator location, 
the tilde represents a skew- symmetrized matrix, and f i is the highest derivative 
of the applied force. Each force derivative resides in a partition of the 
control vector, and is therefore related to the control vector by the selection 
matrix, S L : 

fi = S L u 

On appending the decoupling constraint to the performance index, the resulting 
necessary condition for control optimality is: 

u 2 = -Jr 1 !* 1 * + N T SX * Wk) (33) 
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where 


i 

and k is a 3x1 vector of Lagrange multipliers for the decoupling constraint. 
Substitution of this control expression into the decoupling constraint yields the 
solution for k : 


#c = T X*N r SX) < 34 ) 

Assuming the feedback/feedforward form of Eq. (27), the steady-state gain 
expressions become: 

P = 0 = - PA - A T P + [PS + tf]P' 1 [N I + B T F] - Q ( 35 > 


? = 0 = - [A-BP" 1 (S T ?+N t ) ] t £ - Pd 

where 

R = P + lw J (WR'W)^ 

2 

* — T 
A = A - BR N t 

Q = Q - NR N J 

The final control vector is then given by the expression 

u 2 - 

+ P' 1 ! T [A T -P„SP‘ 1 B T ]' 1 P ss d 


(36) 


(37) 


Extension to Class IV Missions 

The control approach that has been applied to the Class III mission can 
easily be adapted for the Class IV mission. Modifications to the approach 
include the use of multibody dynamics equations instead of Euler's equations for 
the calculation of nominal controls. Partial differentiation of the multibody 
equations would then provide the appropriate costate differential equations. An 
Order (n) algorithm (where n is the number of bodies) would reduce the CPU 
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requirements when the number of bodies is large (Chun, Turner, Frisch, 1989). 
Since the pointing requirements for the appendages would normally be in inertial 
coordinates, and the multibody equations are generally in relative coordinates, 
it is necessary to perform a transformation from inertial to relative coordinates 
in order to specify the desired terminal conditions. 

Another modification to the algorithm would involve expanding the 
disturbance rejection to account for the effects of interaction forces/torques 
as well as the nominal maneuver controls. The perturbation controller would feed 
back the relative appendage motions, in addition to the absolute base body motion 
and the elastic degrees-of-freedom. 
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SIMULATION RESULTS AND OBSERVATIONS 


The EPS large angle slew maneuver is simulated using DISCOS (Dynamics 
Interaction Simulation of Controls and Structure). In this study, the example 
maneuver involves a re-orientation from a frame initially aligned with the 
inertial frame to a frame that is rotated 1 radian about all three Euler axes. 
The Euler sequence is arbitrarily chosen to be 1-2-3. Both torquers-only and 
torquers-plus-force-actuators control are simulated for performance comparisons. 

Slew maneuver times of 6 s and 10 s are considered, which are respectively 
1.45 and 2.42 times the period of the first elastic mode. The purpose of the 
shorter 6 s slew is to investigate the effects on response when there is more 
interaction between the elastic modes and the control torques. 

To evaluate the performance of the different simulation cases and to serve 
as a guideline for the selection of weighting matrices, the following performance 
goals are established. 

(1) Angular errors less than 3% of the desired maneuver angle (1 rad for 
each Euler axis) after 125% of the nominal rigid body maneuver time 
has elapsed. 

(2) Peak elastic displacement at antenna rim less than one tenth of the 
open- loop values after 125% of the nominal rigid body maneuver time 
has elapsed. 

(3) Elastic displacements and control efforts minimized during the 
maneuver. 

These performance goals are based on results of preliminary open- and 
closed-loop simulations. 

The 10 s maneuver results are presented first, followed by the more 
stressing 6 s maneuvers. Open- and closed- loop results are shown for both 
maneuver times. All of the simulated perturbation controllers use a 'constant' 
disturbance rejection, unless otherwise stated. See the section on Disturbance 
Vector on page 18 for details. 

Rigid Body Response to 10 s Open-Loop Torque Profile 

Figure 3 shows the response of a rigid EPS model subjected to the nominal 
rigid body torque profile for a 10 s maneuver. The response should match the 
desired final conditions exactly because the controlled plant is identical to the 
plant model used in generating the nominal torque profile. Indeed, the total 
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angular momentum, Euler angles, angular velocities, and applied torques shown in 
Fig. 3 reveal that the desired final conditions are attained, thus confirming 
that the nominal rigid body torque profile is properly computed. The total 
angular momentum plot shows that the structure's angular momentum increases to 
a maximum and then decreases back to zero at the end of the maneuver as expected 
for a rest-to-rest maneuver. The Euler angle plots show that the structure has 
slewed to 1 radian about all 3 Euler axes at 10 s and stays at that orientation 
throughout the rest of the simulation. The angular velocities are zero at 10 s 
and remain at zero. The applied torques are smooth and continuous because of the 
frequency- shaping formulation, which also allows the torques and their slopes to 
be specified at the initial and final times. These quantities are set to zero for 
continuity and smoothness. The peak torque is about 12,000 N-m about the body 
z axis . 

Flexible Body Response to 10 s Open-Loop Torque Profile 

The first 10 elastic modes are retained in order to observe the response 
of the flexible structure due to application of the rigid body open- loop control 
profile. The results of this case are used to establish a baseline for 
evaluation of the perturbation feedback controller. The Euler angles, body 
angular velocities, and modal amplitudes are shown in Fig. 4. The modal 
amplitude plots show that only the first few modes exhibit any significant amount 
of residual vibration. 

The Euler angle time histories show that there is rigid body drift after 
the end of the maneuver. This is also indicated by the non-zero residuals in all 
components of the angular velocity vector. This rigid body drift is always 
expected when an open-loop control profile, designed with a rigid plant model, 
is applied to a non-rigid plant. 

The elastic deformations at nodes 66 and 69 of the large antenna, and nodes 
85 and 88 of the small antenna are presented in Fig. 5. The elastic deformations 
on the small antenna are seen to be generally a factor of 4-5 smaller than the 
large antenna. All of the responses are dominated by the first and second modes. 

Perturbation Control and Reduced Order Model 

Since the open-loop controller is designed using rigid body dynamics, the 
perturbation feedback control is needed to actively suppress the elastic modes 
and drive the perturbed trajectory as close to the nominal as possible. The 
weight matrix selection process is iterative, and the submatrices, Q and D, are 
constrained to be diagonal so as to limit the number of values to be selected. 
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Figure 5. Elastic Response to 10 



Open-Loop Torque, Physical Deformations 









For the 10 second slew cases when 45 modes are simulated, the results are 
similar to the cases where only the first 10 flexible modes are retained. The 
reduced order perturbation controller is designed using a 10 mode model. 

Perturbation Control - Torauers-Onlv. 10 s Maneuver 

The torquers, which are located at the CG of the vehicle, are primarily 
used to produce the open- loop torques which slew the structure. These torquers 
are also used by the perturbation controller to control the elastic degrees-of- 
freedom as well as deviations from the nominal angle trajectory. This case 
presents the closed- loop response of the flexible plant under perturbation 
feedback control using torquers as the only actuator device. The open- loop 
torque profile is designed for a 10 s slew. For this case, the weighting 
matrices, Q and D, associated with the plant states and control states, 
respectively, are: 
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where I n denotes an nxn identity matrix. The break frequency, Wj, is chosen to 
be 2.5 rad/s. The ordering of the plant states penalized by the weighting matrix 
Q is given by the angular velocities, modal amplitudes and amplitude rates, and 
Euler angles. The translational states are ignored. 

Figure 6 shows the resulting Euler angle histories and torque 
perturbations. The Euler angle histories show that the desired orientation is 
attained after some n damped oscillatory response. There is a 10% overshoot in the 
second Euler angle. The Euler angle errors at 12.5 second are within the 30 mrad 
performance goal specification. The applied torques are computed as the sum of 
the nominal torques (Fig. 3d) plus the torque perturbations (Fig. 6b). The peak 
torque perturbations about each axis are 15-40% of the peak values of the nominal 
open-loop torques. The largest peak torque perturbation is about 3000 N-m about 
the y axis. The elastic deformations at nodes 66, 69, 85 and 88 on the antennas 
are shown in Fig. 7. The elastic deformation of the antennas meets the 
performance goals after the nominal slew maneuver is completed due to the action 
of the perturbation feedback controller. During the slew period, the peak 
elastic deformations at different nodes of the antennas decrease by a factor of 
1.5 to 2 when compared with the open- loop values. 
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Figure 7. Feedback Control for 10 s Slew, Using Torquers Only 






Perturbation Control - Torauers + Force Actuators. 10 s Maneuver 


Since the elastic deformations at the outer rims of the antennas are large, 
it is reasonable to consider placing force actuators at the outer rim to enhance 
to performance of the perturbation feedback controller. The weighting matrices, 
Q and D, associated with the plant states and control states, respectively, are 
as follows: 
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Q a - DIAG [ 25 , 100 ,25,7*7.5,25, 100 ,25,7*7.5] 

In the above expression for Q a , 7*7.5 indicates that 7.5 is repeated 7 times. 
The ordering of the plant states is given by the angular velocities, body 
reference point velocity, modal amplitude and amplitude rates, Euler angles and 
translational degrees of freedom. The diagonal matrix D penalizes the three 
control torques and seven actuator forces. 

The Euler angles, applied forces, applied torques, and torque perturbations 
are presented in Fig. 8. The Euler angle histories show that there is more 
overshoot when compared with the closed-loop torquers-only response (Fig. 6). 
Nevertheless, the Euler angle errors are within the performance goal. The 
actuator forces are also quite large as shown in Fig. 8b, where the peak actuator 
force from one of the actuators is about 70 N. The applied torques and their 
corresponding torque perturbations are larger than that of the torquers-only case 
(see Fig. 6). The peak torque perturbation is about 9,000 N-m. This is an 
indication that there is some degree of interaction between the torquers and the 
force actuators; the overall effect is a further reduction in the elastic 
deformations at the antennas' outer rims. The elastic deformations of nodes 66, 
69, 85, and 89 are shown in Fig. 9. Compared with the open- loop case, the 
response of the large antenna is reduced a factor of 5 and response of the small 
antenna by a factor of 3. 

Disturbance Rejection Feedforward 

Disturbance rejection feedforward is implemented to counteract the elastic 
response induced by the rigid body control law. For comparison, the elastic 
deformations at node 69 of the large antenna are presented in Fig. 10, where the 
disturbance is modeled as: a) no disturbance, b) piecewise constant, and c) time- 
varying. The controller configuration is torquers-only, and the open-loop slew 
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Figure 10. Disturbance Rejection Comparisons 





is 10 s. The weight matrices for the states and controls are given by: 
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It can be seen that the reduction in elastic deformation at this node using 
time-varying disturbance rejection feedforward can be as large as 30% when 
compared with the same case without disturbance rejection. Figure 11 shows the 
Euler angle histories for the case with time-varying disturbance rejection. The 
Euler angle errors of 5 mrad at 12.5 s are well within the 30 mrad performance 
goal. The time-varying disturbance rejection case represents the best results 
obtained for torque-only control. 

Flexible Body Response to 6 s Open-Loop Torque Profile 

The corresponding maneuver frequency is 0.17 Hz, which is about 70% of the 
first natural frequency. The 6 s open-loop torque profile is applied to a 
flexible plant model with all 45 modes retained. The Euler angles and applied 
torque profile are presented in Fig. 12. The first five modal amplitudes and 
elastic deformations at nodes 69 and 88 are presented in Fig. 13. The largest 
slope of the Euler angle errors for this 6 second run is 3-4 times larger than 
the corresponding 10 s slew case (Fig. 4). The peak modal amplitudes are also 
about 3-4 times larger than those for the 10 s slew case, while peak elastic 
deformations at nodes 69 and 88 are 3-5 times larger. 
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Figure 11. Euler Angles for Time-Varying Disturbance Case 


u> 

Ln 



TIME IN SECONOS 


| i>i i | i i i i | i i i i | i i i i ) 

0 5 10 15 20 

R 

TIME IN SECONOS 


Figure 12. Elastic Response to 6 s Open-Loop Torque 
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Perturbation Control - Torauers Only. 6 s Maneuver 


The perturbation controller is simulated with torques only. Linear 
translations are ignored for this case. The weighting matrices used are as 
follows : 
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The break frequency is chosen to be 5 rad/s. 


The simulation results are shown in Figs. 14 and 15. The Euler angles show 
a sluggish behavior. At 7.5 s (125% of the elapsed maneuver time), the Euler 
angle error is about 10-20%, which is outside the performance goal. The peak 
torque perturbation is about 40% of the largest nominal torque. The deformation 
magnitudes at the antenna nodes are about 2-3 times larger than the 10 s 
maneuver, as expected. 


Perturbation Control - Torauers + Force Actuators. 6 s Maneuver 

The closed- loop response of the perturbation feedback with torquers plus 
force actuators are shown in Figs. 16 and 17. The break frequency, Wj, is chosen 

to be 5.0 rad/s. 

Figure 16 indicates that the largest Euler angle error at 7.5 s (125% of 
the maneuver time) is about 7% of 1 rad, which is outside the desired performance 
goal, but better than the case with torquers only. There are also overshoots 
which peak at about t - 5 s. The response is sluggish when compared with a 
similar 10 s slew. A peak force of about 180 N is needed at nodes 60 and 66 
versus 80 N needed at nodes 63 and 69. Modal amplitudes of the first 10 modes 
are presented in Fig. 17. The modal amplitudes of the first 10 modes are 
considerably reduced after the completion of the nominal slew. The peak elastic 
deformation during the slew maneuver is observed to be reduced by a factor of 
about 100 : 1 . 


Summary of Simulation Results 

The perturbation feedback controller provides a stable and satisfactory 
closed- loop slew response. In general, the additional use of force actuators 
yields a slightly better response. However, because of the difficulty in 
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Figure 16. Feedback Control, 6 s Slew, Torquers + Force Actuators 


















selecting the weight matrices to produce comparable performance in both angular 
error and elastic deformation, the difference between torque-only control and 
torque-plus-force control is difficult to quantify. For the range of weighting 
matrices considered, it appears that the benefits of the additional force 
actuators are not significant. 

In general, peak elastic deformations during the slew maneuver can be 
reduced by a factor of 2 when compared with the open- loop response. The typical 
reduction of the elastic deformation during the slew is about a factor of 4. The 
peak actuator force is about 70 N which could make the implementation of force 
actuators on the antennas rim difficult, if not impossible. However, one must 
keep in mind that the EPS structure is not designed for large rapid slews. 
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SLEW CONCLUSIONS AND RECOMMENDATIONS 


Slew Conclusions 


A DISCOS-based multibody simulation has been developed for the EPS model. 
The three -axis frequency- shaped slewing methodology of Chun, Turner, and Juang, 
1987 has been implemented on this model. A baseline set of actuators has been 
defined for EPS involving torquers for the open- loop control, and a combination 
of torquers and force actuators for the perturbation feedback. Both 6 s and 10 
s slew maneuvers have been designed for the EPS. 

Results of the simulation test cases show that the combination of open-loop 
nominal profile, perturbation feedback, and disturbance feedforward performs 
well. The orientation errors and elastic deformations are reduced within a 
reasonable time after the end of the nominal maneuver time. All 10 s maneuver 
simulations meet the performance goal of 3% angular error and 10:1 reduction of 
open-loop deformations. The 6 s example maneuver requires a longer amount of 
time before the angular errors are reduced to that level; on the other hand, the 
elastic deformations are reduced by a factor of 100:1, which is an order of 
magnitude better than the design goal. An adjustment of the weight matrices 
could be made to improve the angular performance while keeping the elastic 
response within the 10:1 reduction goal. 

The simulation results point out the difficulty of selecting proper weight 
matrices for achieving the desired response, particularly when there is a large 
number of states. For simplicity, the weight matrix selection has been limited 
to diagonal elements, even though is known that selection of non-zero off- 
diagonal elements allows greater mobility of the closed-loop poles. The use of 
off-diagonal weighting elements has not been considered for lack of time and 
appropriate software. It is anticipated that with an automated weight selection 
algorithm, all of the closed-loop responses can be tailored to the desired 
performance goals . 

Slew Recommendations 


In this study, the open- loop torques are provided by torque -generating 
devices at the CG of the structure. Future studies could consider the use of 
force actuators to provide part or all of the required torque. With suitable 
allocation of the force levels, this would reduce the amount of excitation to the 
elastic modes of the structure. 

Actuator placement algorithms that account for the effects of the open- loop 
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torques should be used. Improvements to the perturbation feedback controller 
could include the constraining of perturbation forces so that they apply zero net 
torque to the structure and the development of a weight matrix selection 
algorithm to assist in controller design. 

Robustness issues should be examined. These would include robustness to 
modeling errors in the structure and the actuators, and external noises. 

A differential dynamic programming approach could be considered for 
producing a rigorously optimal solution to the open-loop finite-time rigid/flex 
maneuver problem as well as providing the optimal time -varying feedback gains for 
the perturbation controller. Time-varying gains could be fitted to Errors 
coefficients or other basis functions for use in gain scheduling. 

This study demonstrates that the basic idea of open-loop nominal plus 
perturbation feedback is feasible. The choice of an LQR type of perturbation 
feedback is but one of many possible controller approaches. Future studies 
should consider different types of perturbation feedback controllers, depending 
on the mission requirements and computational capabilities. 


INTEGRATED DESIGN: TECHNICAL AND SOFTWARE DEVELOPMENT ISSUES 


Integrated design seeks to bring together several technical disciplines in 
a common work environment to model, analyze, predict performance, and optimize 
the controlled behavior of large flexible space systems. The analysis problem 
is multidisciplinary in nature. A complete investigation of the expected on- 
orbit performance of a mission concept can require a structural model, thermal 
model, disturbance model(s), control/estimation model , performance metric models, 
optical and/or radiometric models, dynamics models, ground test models and 
potentially other models. This process starts during the preliminary design 
phase for a new concept and continues throughout the life of a mission. 
Traditional spacecraft design approaches have developed approaches that have 
segregated activities along discipline lines. With the next generation of 
planned satellites there is an increasing need to closely coordinate the 
activities of different design groups in order to insure that the performance 
goals for the next generation of large lightweight joint-dominated systems can 
be achieved while reducing the mission risk to unanticipated design problems. 

Design Issues for Large Space Structures 

The design of Large Space Structures(LSS) is a multidisciplinary activity 
which embraces structures, structural dynamics, control, and multibody dynamics. 
The process starts with the identification of a new mission and the performance 
requirements required for achieving the science and engineering objectives. As 
science missions become more ambitious the spacecraft design community is forced 
to develop large systems in order to achieve the sensitivity and resolution 
requirements for analyzing the mission observables. An important byproduct of 
the need for large systems is that launch mass limits force the structural design 
towards lightweight low frequency structures. The performance tolerances tighten 
because larger systems have inherently greater resolving power. The resolving 
power of the instruments places constraints on the allowed levels of jitter 
motion for the structure. The characterization of the coupling effects between 
vehicle size and instrument performance is complicated because the effects are 
functionally dependent on the operational wavelength of the science instrument, 
the inherent resolving power of the device, the physical dimensions of the 
operational platform, and the signal processing capabilities available for 
analyzing the data. Large low-frequency systems can be shown to be highly 
susceptible to experiencing an induced response, and that the susceptibility 
increases nonlinear ly with increasing vehicle size. Fortunately, heuristic 
criteria can be established that can be used as guidelines for determining when 
a new structural concept may require active control means for achieving the 
mission objectives. When an application has been identified where control is 
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required, then advanced technology can be beneficially applied and may in fact 
represent an enabling capability. 

Genesis of the Interaction Problem 

Mathematically, the need for active control arises when the bandwidths for 
the structural frequencies, the disturbances and the control are either close or 
overlapping (see Figure 18). When these bandwidths are close or overlapping a 
real possibility exists for exciting a resonant structural response which can 
either degrade the performance or damage the vehicle. The complexity of the 
control problem is assessed by evaluating the open- loop performance of the 
structure and determining whether any induced responses exceed the structural 
design performance specifications. When the response exceeds the design 
specifications then some form of control must be invoked in order to achieve 
acceptable system- level performance. There are four levels of control design 
complexity which characterize the difficulty of the problem, as listed in 
Table 1. 

Table 1: CONTROL SYSTEM DESIGN COMPLEXITY 


Problem 

Type 

Performance 

Spec.* 

(Range) 

Control Benefit 

I 

lO^lO 1 

Enhancing Capability 

II 

lO^lO 4 

Enabling Capability 

III 

10-10 6 

Enabling Capability/ 
Multiple Technologies 
Required 

IV 

10>10 6 

Beyond the State -of- the 


* (Open-Loop Performance)/(Performance Goals) 


Type I missions represent enhancing capabilities because many potential 
solution strategies are available, including: slight re-designs of the 
structure, added damping treatments, and system/component level isolation 
concepts . 

Type II missions represent enabling capabilities because both classical and 
modern control approaches can be used to achieve the required levels of the 
system performance. As discussed in Turner, 1990, ACOSS-type control can be used 
for broadband disturbance inputs that induce open- loop responses exceeding the 
performance specifications in the range of [10 l , 10 2 ] . This capability is 
consistent with the goal of many ground-based experimental efforts, though to 
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Figure 18. CSI Problem Complexity 
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date no results have been reported that corroborate the theoretical predictions 
throughout the performance range of interest. For response ratios in the range 
of [10 2 , 10 4 ] there are theoretical predictions which indicate that single- 
frequency disturbances can be suppressed using ACOSS technology. 

Type III missions are enabling but require hierarchical approaches which 
combine passive damping treatments, active isolation concepts, broadband/single 
frequency modern control approaches, and optimization methodologies for obtaining 
a workable design. 

Type IV missions are characterized by control requirements which exceed the 
state-of-the-art in one or more critical design technology areas. 

Types II and III represent the applications where advanced design 
capabilities can be beneficially used for solving complex structures and control 
design problems. The technology issues which impact the development of an 
integrated analysis capability are presented in the following sections. 

Control/Structure Interaction Technology 

Control/Structure Interaction (CSI) technology provides a means for 
suppressing undesirable structural behaviors and achieving the mission objectives 
by either passive or active means. CSI is generally an enabling capability for 
systems which are either very large or subject to tight internal alignment 
constraints typical of optical systems. 

In traditional spacecraft design approaches the structures and control 
teams have minimal levels of interaction during the design effort (see Figure 
20). This approach works fine as long as the structure is basically rigid and 
the attitude and subsystem pointing controls and disturbances cannot induce a 
structural response. A key assumption of this approach is that small changes in 
either structural or control design will overcome any potential problems which 
may arise. As vehicles become larger, however, the need for coordination between 
design teams increases because many potential problems can be easily fixed during 
the early phases of a project. On the other hand, if the problems are only 
uncovered late in a project, the impact on time and cost can be significant. 
Moreover, failure to detect CSI related design flaws can lead to either greatly 
degraded system performance or catastrophic system failure. 

A fully integrated analysis approach can offer a number of distinct 
advantages for a spacecraft design group. The key objectives and payoffs which 
result form an integrated software design effort. 
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Technical Issues 

Potential Concerns ' 

Resolution/Tradeoff 

Modeling Accuracy 

•Nonlinear Effects 
•Fluid Motion 
•Stability 
• Performance 

•Material Nonlinear Joint Fluid 
Motion Models Increased Fidelity 
vs. Computational Cost 
•Ground Testing 
•System ID, Adaptive Control 

Control Law Design 

•ROM Spillover 
•Robustness 
• Performance 

•Improved Model Selection/ROM 
•Robust Design Techniques 
•Nonconservative Analysis 
Techniques 

•Distributed Control Approaches 
•Adaptive Control 

Sensors & Actuators 

•Poor Performance Due to 
Bandwidth, Saturation, 
Resolution, Noise, Nonlinearity, 
Placement 
•Failures 

•New Technology (Fiber Optics, 
Electro-Optics, Piezoelectrics) 
•Optimal Placement Methods 
•Testing and Refinement 
•Redundancy Management 

System Identification 

•Resolution of Closely Spaced 
Modes 

•Computational Complexity 
•Sensor Resolution and Bandwidth 
•Nonlinearities 

•Stochastic Methods 
•Adaptive Methods 
•Parallel Processing 

System Validation 

•Testing in 1-G 
•Facilities for LSS 
•Suspension & Scaling 
•Cost 

•Develop Space-Validated 
Analytical Tools 
•Novel, Active 0-G Suspension 
•Scaling Laws 
•New Facilities 


Figure 19. Key CSI Technical Issues 
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Obi ectives 

• Automated tool for preliminary and advanced structural 

and control design 

• Emphasis on multidisciplinary aspects of the design task 

Payoffs 

• Improved design team communication 

• Enhanced designs 

• Testbed for evaluating emerging analytical design 

methodologies 

• Reduce CSI risks 

The structural behavior is measured relative to the mission performance 
goals. To measure the system performance a model is developed which consists of 
a normal mode model (e.g., typically obtained from a finite element model ) and 
a disturbance model which consists of both on-board and environmental mechanical 
loads. By simulating the performance of the system without CSI control 
technology the complexity of the problem can be assessed. For example, if the 
line-of-sight pointing error and jitter exceed the allowed tolerances for 
achieving the science and engineering goals, several options are available for 
carrying out corrective action, including: 

• Structural Redesign 

• Advanced Materials/Damping Treatments 

• System/Component Isolation Concepts 

• Vibration Control 

-Shape Control 
-Figure Control 
-Slew Control 
-Vibration Suppression 

An integrated design environment provides an opportunity for evaluating new 
and emerging control design approaches in a standardized evaluation context. 
This approach is ideal for minimizing the time required for transferring new 
technology development efforts from the research community to the spacecraft 
design community. 

The key technical issues which affect the complexity of CSI design efforts 
are listed as follows: 

Modeling Accuracy 
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Control Law Design 
Sensors and Actuators 
System Identification 
System Validation 

Each of these issues is complicated because of explicit and implicit coupling 
effects that arise through modeling idealizations and approximations, algorithmic 
limitations, as well as hardware limitations and manufacturing errors (see Figure 
19). 


CSI Technical Issues: Modeling Accuracy 

Modeling accuracy defines how well the physical parameters describing the 
structural model are known. These parameters can include: geometry, mass 

distribution, stiffness distribution, frequencies, mode shapes, joint behavior, 
control system dynamics, and disturbance dynamics. This issue is important 
because errors in the mathematical model directly impact the stability, 
robustness, and performance characteristics of the control design effort. When 
modern multi-input and multi-output control approaches are invoked it is well- 
known that high levels of system performance can only be maintained when both the 
model for the structure and the estimated structural behaviors are accurate. An 
early and very important step in the model development process consists of 
selecting a suitable reduced-order model from the finite element description of 
the structural behavior. This step is important because it involves trade-offs 
between model fidelity, computational complexity, system performance and 
stability. Several of the currently popular reduced-order model selection 
methodologies are presented as follows: 

Gain Factor Analysis 
Singular Value Decomposition 
Internal Balancing 
Modal Cost Analysis 
Chained Aggregation 

A discussion of these techniques can be found in Turner, 1990. 

CSI Technical Issues: Control Law Design 

A control law design is basically a mathematical process whereby a strategy 
is developed for using estimates for the instantaneous system state to correct 
errors in the system performance. Either classical or modern control approaches 
can be used. However, only the modern control approaches have the theoretical 
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Figure 20. Conventional System Design Flow Diagram 
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capability for handling multiple inputs and outputs simultaneously. A partial 
listing of the currently available control design methodologies are presented in 
what follows: 

Direct Output Feedback 
Parameter Optimization 
Orthogonal Filter 
Positivity/Characteristic Gain 
High Authority/Low Authority Control 
Model Error Sensitivity Suppression 
LQG/LTR 
H - , H 2 , L" 

Maximum Entropy 
Adaptive Methods 

These methods are discussed in Turner, 1990. 

CSI Technical Issues: Sensors and Actuators 


Sensors are used to obtain measurements of the system state that can 
consist of translational and rotational position and velocity data. The 
actuators provide the means for carrying out force and torque commands for 
achieving closed-loop control objectives. The limiting factors which affect the 
accuracy of these devices are bandwidth, resolution, biases and scale factors for 
the inputs and outputs. These components act as the ultimate constraints on the 
achievable performance for any system. 

CSI Technical Issues: System Identification 

Because any control system design is based on various modeling assumptions 
a high performance control system can only be expected to operate successfully 
when the model parameters are well-known. The approach used for determining the 
model parameters form observations of the system response is system 
identification. Through analysis, system identification refines parameter 
estimates so as to make the model predictions more closely agree with the 
observed behavior. This technology becomes particularly important as the 
performance specifications tighten the limits on the allowed motions, thereby 
forcing the control system to perform extremely accurate feedback commands. 
Several of the currently available algorithms for carrying out the system 
identification function are listed as follows (Turner, 1990): 
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Cross -Correlation 
Transfer Functions 
Kalman- Filtering 
Lattice Filter 
ERA 

Maximum Likelihood Method 
ITD 


CSI Technical Issues: System Validation 

Both ground- and space -based system validations are required in order to 
ensure that the performance specifications are achieved and that the system is 
stable when operating in the presence of realistic disturbance sources. Tests 
of the actual operational hardware provides the only reliable methods for 
assessing the behavior of complex systems. The tests also provide ideal 
opportunities for validating analysis tools and improving one's understanding of 
the modeling process. When system structural concepts become too large to 
experimentally validate on the ground, one is forced to depend exclusively on 
analysis tools. For this class of applications there is a clear need for 
developing, testing, and validating software analysis tools. 
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INTEGRATED DESIGN APPROACH 


Integrated design implies that the needs of several technology areas are 
considered simultaneously. For example, when the needs of structures and control 
are considered, one can seek a minimum weight structure where the control/ 
estimation/system identification gains are determined subject to constraints on 
allowed subsystem motions, closed- loop pole locations, and numbers of sensors and 
actuators. For an integrated design approach to be effective the desired 
engineering approach consists of having all computational tasks carried out in 
a common simulation environment (see Figure 21) . This approach ensures the 
consistency of data links from the structures design, control design, thermal 
design, optics/radiometric design, multibody models , optimization , and system- 
level evaluations. An integrated design strategy seeks to exploit a team 
approach for addressing design problems, thereby minimizing the traditional 
institutional barriers which act to segregate design group members along 
discipline lines. 

Review of Research in Integrated Design 

Integrated design during the last eight years has been used to refer to 
simultaneous or dual optimization studies where the design variables have 
consisted of both structural and control parameters. A major goal of an 
integrated design effort is to determine a better design for vehicles, where the 
success of a design is measured by evaluating performance measures such as 
minimum time, minimum fuel, minimum control effort, minimum weight, and minimum 
cost. 


Mathematical optimization approaches are used for successively refining the 
structure and control designs, because of the multidimensional nature of the 
problem. A generic description of the mathematical statement of an optimization 
problem is presented in the following section. 

Generalized Kuhn-Tucker Conditions 

An optimization problem is defined by specifying the free variables which 
can be adjusted in a successive approximation algorithm that seeks to improve 
some measure of the system performance. To ensure that a physically meaningful 
set of design variables results from the optimization process one must specify 
constraints. For example, cross sectional beam areas must always be positive 
real numbers. The searching algorithms are provided measures of the system 
performance, such a quadratic performance indices which penalize structural 
response and control effort. By measuring the sensitivity of the performance 
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measures, one can define incremental changes in the design variables that 
minimize (maximize) the performance measures locally. Mathematically the 
optimization process can be described as follows: 

Find D t , i - 1 n, such that (D A a Design Variable) 

Minimize (Maximize) F(D 1 ,...,D n ) (Scalar or Vector) 


Subject to the Constraints 

hj(D 1 , . . . , D n ) - 0, j - 1 p (Equality Constraints) 

gk( D i« • • • » D n ) - 0. k - 1 q (Inequality Constraints) 


Form the Augmented Objective Function 




A' S' 

F(D 1 , , D a ) > • • • < t <7* > • • • » D a ) + 5’jc) 

jc-1 


(38) 


where 

Aj Denotes the Equality Constraint Lagrange Multiplier 
Denotes the Inequality Constraint Lagrange Multiplier 
S k 2 Denotes the Inequality Constraint Slack Variable 

The necessary conditions which define the mathematical requirements for 
determining an optimal solution are: 

P 7 

dF/dD^V Xjdhj/dDi + Y, , (i*l* . . . , n) 

jZ\ *•! 


where 


/ijj - 0 if g k (*) > 0 and p* > 0 if g k - 0 

The basic mathematical problem statement is highly nonlinear and successive 
approximation strategies provide the only rational approach for systematically 
refining estimates for the optimal set of design variables. Numerical 
optimization problems are difficult for two reasons. First, initial starting 
guesses must be specified which must satisfy the constraint conditions for the 
algorithm to start with a feasible solution. Second, the analyst must be 
constantly alert to the possibility that multiple minima (maxima) can attract the 
optimization algorithms, leading to converged solutions that may not be globally 
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optimal . 


Optimization based design techniques have been developed for many years 
(Crane, Hillstrom and Minkoff, 1980; Arora, Thanedar and Tseng, 1985; Grandhi, 
Thareja and Haftka, 1985; Gill, Hammarling, Murray, and Wright, 1986; Gill, 
Murray, Saunders and Wright, 1986; Gill, Murray and Wright, 1981; Polak, Siegel, 
Wuu, Nye and Mayne, 1985; Palacios -Gomez , 1980; Vanderplaats , 1973; Goldfarb, 
1969; Venkayya and Tischler, 1978, 1979). At this time there are several 
commercially available tools which have implemented several of the 
multidimensional search algorithms. They include: 

Recursive Quadratic Programming 
Homotopy- Based Optimization Methods 
Gradient Projection 
Reduced Gradient Methods 
Feasible Directions Methods 
Optimality Criteria Methods 
Sequential Linear Programming Methods 
Projection Methods 

Sequential Unconstrained Minimization Techniques 
Multiplier Methods (Augmented Lagrangian) 

The usefulness of any of these methods depends on the special features of the 
problem at hand. Some of the methods make exclusive use of first-order 
sensitivity data, others either compute exact expressions for the second-order 
sensitivity data or approximate this information by various algorithmic 
strategies (Heylen, 1986; Rudisill, 1974; Fox and Kapoor, 1968; Plaut and 
Huseyin, 1973; Haug and Rousselet, 1980a, b). Common to each of the methods is 
the need for defining measures or performance indices for measuring the 
effectiveness of parameter refinement strategies, evaluating any constraints 
which limit the allowed range of acceptable solutions, computing the sensitivity 
of the performance measures and the active constraint sets to variations in the 
current estimates for the design variables. 

The calculation of the sensitivity data represents a computationally 
intensive part of any optimization-based study. There are two frequently used 
methods for generating the required partials: finite difference techniques and 

analytic partials. The finite difference technique is conceptually the simplest 
approach because one uses numerical perturbation techniques to generate 
approximations for the required partial derivatives (Rudisill, 1974; Cronin, 
1988; Mills-Curran, Lust and Schmit, 1983; Lust and Schmit, 1986; Fox and Kapoor, 
1968; Prasad, 1984; Fleury and Schmit, 1980; Palacios -Gomez, 1980; Vanderplaats, 
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1984; Goldfarb, 1969; Plaut and Hus ey in, 1973). A potential disadvantage of the 
finite difference technique is that the numerical procedure can be iterative and 
expensive if the model itself is computationally intensive to evaluate. The 
analytic partials approach assumes that closed- form expressions can be developed 
which can be used to provide the sensitivity to machine precision. When 
practical the analytic partials approach is preferred because all the data can 
be generated with a singe operation. 

Optimization for Structural Applications 

There is an extensive literature which describes application of 
optimization-based methodologies for structural applications (Fraser, 1985; 
Grandhi, Thareja and Haftka, 1985; Choi, Haug and Seong, 1983; Zeischsha, 
Leuridan, Storrer and Vandeurzen, 1988; Heylen, 1986; Cronin, 1988; Hou and Twu, 
1986; Woo, 1987; Khan, Thornton and Willmert, 1982; Yoshida and Vanderplaats , 
1985; Rubin, 1970; Fleury and Schmit, 1980; Pedersen, 1982; Shiau and Chang, 
1988; Kamat, Venkayya and Khot, 1983; Khot, 1985; Grandhi and Venkayya, 1987; 
Shin, Plaut and Haftka, 1987; Thareja and Haftka, 1986; Canfield, Grandhi and 
Venkayya, 1986; Venkayya and Tischler, 1978). Typical design approaches define 
performance measure(s) and constraints (geometrical-manufacturing limits, 
practicality, behavioral -pole placement, peak deflections) . Typical structural 
performance indices include: 

Cost 

Weight 

Stress 

Stiffness 

Deflection 

Frequency 

Mode Shape 

Fatigue Life 

Reliability 

Many times it may be necessary to consider several performance indices 
simultaneously in order to obtain a practical design. A standard approach is to 
consider a linear combination of performance indices that are multiplied by 
weighting coefficients which are assigned based on some heuristic criteria. This 
approach is capable of producing numerical results that satisfy the necessary 
conditions for the optimal solution. Nevertheless, there are many tradeoffs 
between the different objectives that can not be fully investigated via this 
approach. The alternative to the scalarization of several performance indices 
is to consider vector performance indices. This approach, however, has received 
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little attention in the technical community. 


Typical structural design variables for an optimization process include: 

• Cross-sectional areas 

• Thickness 

• Inertial & Torsional Properties 

• Topological (Members & Joints) 

-Number 

-Sequence 

-Connectivity 

• Geometrical (Coordinates of Joints) 

-x,y,z ,<i>,9 ,ij> 

• Shape 

-Cross-Sectional 

-Boundary 

• Material 

-Fiber Orientation 
-Passive Damping Distribution 

• Loading 

-Location 

-Number 

-Type of Support 
•External Loads 


Optimization for Control Applications 

In any control scheme there are free parameters which can be adjusted in 
order to achieve desired levels of system performance. Indeed, in modern control 
approaches the necessary conditions for the state and costate differential 
equations require that the elements of several matrices be specified, leading to 
a parameter optimization problem at the very start of the problem. Because the 
necessary conditions for an optimal control solution are naturally described by 
coupled systems of nonlinear differential equations, where the known boundary 
conditions are split between the initial time and final time, many first- and 
second-order optimization methods have been developed for solving these problems. 
The trajectory control problems are typically open- loop control problems, though 
several second-order methods also produce time-varying feedback gains. The use 
of optimization-based methods has been a natural part of the solution process for 
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control problems that have been defined by integral-variational formulations. 
(Boyd, Balakrishnan, Barrat, Khraishi, Li, Meyer and Norman, 1988; Desoer and 
Gustafson, 1984; Doyle, 1982; Doyle and Stein, 1981; Francis and Zames, 1984; Ly, 
Bryson and Cannon, 1983; Polak, Siegel, Wuu, Nye andMayne, 1985; Safonov, 1986; 
Garibotti, 1984; Walsh, 1987; Manning and Schmit, (1987). 

Many performance measures can serve to define optimal control problems. 
One of the great benefits of numerical optimization methodologies is that 
classical design criteria can be combined with modern multi- input/multi -output 
designs in a systematic way. As in the case of structural when scaler 
combinations of different objectives are combined many potential options are 
lost. Alternatively, vector objective functions may be the best approach for 
dealing with conflicting design objectives. Typical control performance goals 
measures and technologies are listed as follows: 

• Pole locations 

-Control Gains 
-Estimator Gains 

• Performance Index(es) 

• Control Type 

-Regulator 

-Tracking 

-Disturbance Accommodation 
-Frequency Shaping 
-Differential Dynamic Programming 

• Controllability 

• Observability 

• Robustness 

• Transfer Function Behavior 

• Loop Stability 

• Frequency- Domain Loop Shaping 

• Design Criteria 

-Stability Margins 
-Bandwidth 
-Damping 
-Overshoot 

The design variables for the control design can include both algorithm 
related parameters and hardware design and placement issues. Typical control 
design variables are listed as follows: 
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Sensor/Actuator Placement 
Weight Matrices 




-State 

-Control 

-Terminal Penalties 
-Estimator 

-Feedback/Feedforward Gains 
Optimization for Integrated Design Applications 

The application of optimization-based techniques for simultaneous structure 
and control design began with the work of Hale, Lisowski and Dahl, 1982, in which 
they considered a single-axis slewing control problem for a flexible structure. 
The basic structure consisted of a central hub with four symmetrically attached 
elastic appendages. They considered the problem of minimizing a quadratic 
performance index. Their control formulation consisted of an open- loop maneuver 
profile. The design variables for the structure consisted of the boom taper and 
the design variable for the control consisted of the state and control weighting 
matrices. The coupling between the control and structural problems arises 
through the closed- form expression which is obtained for the performance index, 
because the performance index is minimized subject to the equations of motion as 
a differential equation constraint. A gradient projection algorithm is used to 
refine the approximations for the design variables. This work stimulated many 
additional investigations into the area of integrated design methodology. A 
brief summary of other related work is presented in Table 2. The authors listed 
in Table 2 represent the individuals who have made significant and innovative 
contributions to the literature. 
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Table 2. RESEARCH IN INTEGRATED DESIGN 


Author (s) 


Problem Description 


Salama, Hamidi and Demsetz, 

1984 

Messac and Turner, 1984 

Haftka, Martinovic and Hallauer, 

1985 

Junkins, Bodden and Turner, 1984 

Miller and Shim, 1986 

Khot, Grandhi and Venkayya, 1986- 
present 

Adamian and Gibson, 1986 

Cooper, Young and Sutter, 1986 
Kyong and Junkins, 1988 

Lust and Schmit, 1988 

Gustafson, Aswani, Doran and Tseng, 
1985 

Miller and Shim, 1986 
Belvin and Park, 1988 & 1990 
Onoda and Hafka, 1987 


Feedback Control , LQG, Second-Order 
Newton-Raphson Iteration, Steady- 
State Riccati Equation 

LQG, Steady-State and Finite-Time 
Control, Single-Axis Slew Control, 
Minimum Weight Structure 

Grillage Control, Steady-State and 
Control, Simplified Structural 
Changes for Control, Theory and 
Experiment 

Generalized Quadratic Performance 
Index, Slewing, Homotopy Minimum 
Norm Optimization 

Truss Structure, Gradient Search, 
Steady-State 

Steady-State Control, Minimum 
Weight, Gain Norm Constraints 

LQG, Minimum Weight, Control 
Robustness, ADS Optimization 

IMAT Based Structure/Control Design 

Open-Loop, LQG, Sensors and 
Actuator Placement, Homotopy/ 
Simplex Optimization 

Direct Output Feedback, CONMIN 

Continuum and Frequency Domain 
Design Techniques 

Minimum Weight and Minimum Control 
Energy, Gradient Optimization 

LQG, ADS Optimization, LASS, 
ORACLS, CS 3 Prototype System 

Minimum Weight and Cost 
Minimization 
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There are two classes of spacecraft operational conditions which require 
different control approaches. First, finite time transient problems which are 
typical of slewing control problems. These problems are not likely to induce a 
resonant response in the structural behavior because the slew disturbance stops 
at the end of the maneuver. Second, infinite time problems which correspond to 
steady-state vehicle operations where a persistent disturbance source(s) acts to 
induce a structural response which must be suppressed by active means. The 
issues of control system stability and robustness are more critical for problems 
of the this type. Most of the research activity listed in Table 2 has been 
focused on the finite time slewing control class of optimization problems. To 
date there have been no realistic infinite time problems and solution approaches 
discussed in the technical literature. This issue is important because the 
majority of real spacecraft applications can be expected to be of infinite time 
variety, hence, a need remains for exploring the computational and algorithmic 
issues associated with this class of technical problems. 

Several general conclusions can be reached regarding the research conducted 
to date on the problem of integrated design. First, this work has demonstrated 
the basic feasibility for conducting fully integrated design efforts. Second, 
much of the effort of the individual researchers has focused on obtaining the 
analytic partial derivatives required for the optimization algorithms. This 
information has been widely disseminated in the literature but no general-purpose 
capabilities exist at any one organization. In problems where researchers have 
used approximate finite difference techniques to compute the sensitivity partial 
derivatives it has been consistently reported that the structural partial 
derivative calculations dominate the computer run times. As a result, this 
observation has stimulated much work in the area of developing analytic partials 
for the structural models, using typical eigenvalue/eigenvector data obtained 
from finite element programs such as NASTRAN. Third, there has been no 
community-wide convergence to a best method for attacking this class of problems, 
though several methods appear to work equally well. Fourth, no large -order 
problems have been exercised to date (i.e., less than 100 degrees of freedom). 
Fifth, some simplified approaches have been presented and appear to be useful for 
limited classes of problems. Sixth, no theoretical obstacles have been 
encountered which can limit the extension of existing results to problems of 
arbitrary complexity. 

Integrated Analysis Software Design Capability 

The development of an integrated analysis capability must first identify 
the expected user group for the tool in order to ensure that the flexibility and 
generality required in an operational spacecraft design environment is properly 
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developed. As presented in Figure 22 there are three levels of user groups which 
can be expected to interface with an integrated design tool. First, system 
engineers whose basic interest is in conducting system- level performance 
assessments and trade studies, particularly during preliminary design studies. 
To support these users the integration software capabilities across all technical 
disciplines is an absolute requirement. Here algorithm robustness is more 
critical than efficiency because these engineers will generally be interested in 
low- dimensioned problems. These users would not be expected to understand the 
technical underpinnings of all of the software packages they must use. As a 
result, the integrated analysis capability must be design to function as a black 
box for these users, because they would never be expected to modify the core 
analysis software in any way. Second, discipline specialists whose primary 
interest would be in conducting detailed trade studies in an individual technical 
area. Here only a medium-level of problem integration is required across the 
supporting technical disciplines. This class of users would be expected to 
conduct stressing performance assessments with high-fidelity/high-dimensioned 
problems in their technical areas of interest. These users would not be expected 
to modify the core analysis software in order to carry out their studies. Third, 
basic researchers whose major goal is the develop and testing of advanced 
capabilities in their technical disciplines. Here only a low-level of problem 
integration is required across the supporting technical disciplines. Typical 
problems can be expected to be of high dimension and representing state-of-the- 
art capabilities. These users would be expected to integrate their software 
developments within the integrated design capability. Clearly the needs of the 
potential user groups are very different but they can all be addressed within a 
common framework. 
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Figure 22. Multi-Level Functionality User Groups for an Integrated Design 
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PROGRAM PLAN FOR DEVELOPING AN INTEGRATED DESIGN TOOL 


Development Options 

There are three development options which should be considered for building 
an integrated design tool. First, one could expand on the multidisciplinary 
optimization efforts that have been under development by either Sobieski, 1982, 
1984; or Venkkaya, 1978, 1979, 1987. This approach would ensure technical 

impeccability of the analysis tools. Unfortunately, these tools have been 

designed for use by researchers in the field of optimization and would therefore 
be difficult for generalists to use. Second, one could seek to build the 
capability from scratch, however. This approach is potentially too costly and 
does not leverage the previous significant government investments in analysis 
capabilities. Third, one can build on existing multidisciplinary tools such as 
IMAT(LaRC), IAC(GSFC), and ISM(AFWL) . IMAT exists at LaRC and provides 
connections for structures and control through a common data base design. IAC 
provides a general architecture for connecting many analysis capabilities but 
represents a precursor to ISM. ISM has the most promise but is currently 
beginning pre-beta testing and will not be available for some time. As a result, 
for the near term the best approach is to build a prototype capability based on 
in-house tools available at LaRC. Over time the supporting capabilities 
developed for the prototype tool could be transferred to ISM for building the 
general-purpose tool. Figure 23 presents a potential structure for a long-term 
integrated design development effort. The basic tool would have an interface 
which accesses ISM. ISM would act as a database manager for controlling the 
internal flow of information between the different technology computational 
modules. The key point to observe is that the optimization, control, structures, 
thermal, optics, radar, and associated sensitivity partial derivatives all share 
common databases, workspace, host files, system utilities, and compiler 
optimization. Through the user interface/executive system, users can control the 
level of integration across the available discipline tools and the level of 
complexity required within individual technology areas. The identification of 
the multiple levels of users needs ensures that a broad base of engineers and 
scientists can fully exploit the unique capabilities that such a tool can 
provide. 

A proposed approach for developing the prototpye capability follows. 
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Prototype Integrated Design Capability 


There are several early design considerations which must be taken into 
account. These include: 

• Modular Design 

-To Readily Accommodate New/Improved Tools 

• Interfaces to Existing Design Tools 

- IAC 
-ISM 
-ASTROS 
-MATRIXx 
-ADS 

• Development of an Application Library 

-Standardized Test Problems & Results 
-Archive Demonstrated SOA Capabilities 

• Planning for Parallel Research & Hardware Validation Efforts 

-Structures 

-Control 

-Optimization 

-Algorithm Development 

-Parallel Computational Architectures 

• Target Computer Platforms 

-High Performance 386 & 486 PC's 

-Workstations 

-Mainframes 

-Parallel Machines 

- Supercomputers 


Beyond these top-level design considerations the specific ingredients of 
a prototype integrated design tool must be defined. To make full use of LaRC's 
existing capabilities, the following approach is proposed for developing a near- 
term in-house prototype capability: 


• Use IMAT as the Core Database Manager 

-MSC/NASTRAN ( Structures Finite Element Code) 
-MATRIXx ( Classical & Modern Control Design Code) 

• Use ADS as the Optimization Tool 

-Industry Standard Software 

• Develop Finite Difference Methodologies for Partials 

-Provides Near-Term Operational Capability 

• Develop Analytical Partial Models for Specified 
Controls and Structures Applications 
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-Implement LQG and Pole Placement Control Algorithms 
-Develop Library of Analytical Partials for Future 
Applications 

• Develop the Prototype Program Executive 

• Validate the Prototype Tool by Analyzing Previously 
Published Applications 

• Develop Ground Test Article for Designing, Testing, and 
Validating the Predictions of the Integrated Design Tool 

The architecture for the proposed system is presented in Figure 24. The 
basic idea is to use the existing in-house capability to the greatest extent 
possible. This approach combines the structural analysis, control synthesis, and 
response predictions that exist through IMAT. The new software to be developed 
will consist of the integrated design interface/executive, ADS, structural 
partials, control partials, and response partials, where the partials are 
calculated initially via finite difference approaches. Under model development 
it is recommended that standard application problems such as the Draper/RPL 
model, planar truss, Class II, and Class IV models be developed. The Class IV 
problems will represent a stressing state-of-the-art application which will 
require that IMAT be modified to run a multibody program such as DISCOS in order 
to generate the platform and robot manipulator trajectory dynamics and loads. 

Tasks and Resources Required to Develop Prototype Tool 

Figure 25 presents a 24-month program for developing an in-house LaRC 
prototype integrated design capability. The first year builds and tests the 
basic connections between software tools and implements several classical and 
modern control approaches. Preliminary integrated design results will be 
available at the end of the first year assuming that the DARPA/RPL model is used. 
A task is also identified for planning the transition from the prototype tool to 
an ISM-based capability in the future. This future capability will be fully 
capable of analyzing and designing any foreseeable system that NASA can envision 
in the next ten to twenty years. Manpower estimates are also provided for 
carrying out the individual tasks. The entire prototype development effort is 
estimated to require approximately four man-years to complete. 
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Figure 24 Prototype Software Development Plan 
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